Enhanced vortex heat conductance in mesoscopic superconductors 
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Electronic heat transport along the flux lines in a long ballistic mesoscopic superconductor cylin- 
der with a radius of the order of several coherence lengths is investigated theoretically using both 
semiclassical approach and the full quantum-mechanical analysis of the Bogoliubov-de Gennes equa- 
tions. The semiclassical approach is constructed analogously to the Landauer transport theory in 
mesoscopic conductors employing the idea that heat is carried by the quasiparticle modes propa- 
gating along the vortex core. We show that the vortex heat conductance in a mesoscopic sample 
is strongly enhanced as compared to its value for a bulk superconductor; it grows as the cylinder 
radius decreases. This unusual behavior results from a strongly increased number of single-particle 
transport modes due to giant mesoscopic oscillations of energy levels, which originate from the in- 
terplay between the Andreev reflection at the vortex core boundary and the normal reflection at 
the sample edge. We derive the exact quantum-mechanical expression for the heat conductance and 
solve the Bogoliubov-de Gennes equations numerically. The results of numerical computations con- 
firm the qualitative Landauer-type picture and allow us to take into account the partial reflections 
of excitations. We analyze the effect of surface imperfections on the spectrum of core excitations. 
We show that the giant oscillations of core levels and thus the essential features of the heat transport 
characteristic to ideal mesoscopic samples hold for a broad class of surface imperfections as well. 

PACS numbers: 74.78.Na, 74.25.Fy, 74.25. Op, 74.45. +c 



I. INTRODUCTION 

Low-temperature measurements of the heat transport 
in the mixed state of type-II superconductors is a unique 
tool providing information about characteristics of the 
quasiparticle excitations and superconducting gap struc- 
ture. When the external magnetic field B is applied, 
thermal transport in the direction of field is dominated 
by the quasiparticle waves propagating along the vor- 
tex cores. Thermal conductivity in the mixed state is 
a longstanding problem considered in many classical pa- 
pers (see, e.g., Refs. 00 where the situation near the 
upper critical field H c2 was discussed). In dirty super- 
conductors, the electron contribution to the thermal con- 
ductance along the vortices is proportional to the area 
occupied by the cores k(B) ~ (B / H C 2)kn , where kjv is 
the electron thermal conductance in the normal statei 3 ^^ 
In clean superconductors this simple estimate fails to 
describe the experimental data-£ the thermal conduc- 
tance along the well-separated vortices in s-wave super- 
conductors appears to be two orders of magnitude smaller 
than predicted. This discrepancy was first explained in 
Ref. Q as resulting from the very small group velocity 
of the Caroli-de Gennes-Matricon (CdGM) states^ lo- 
calized inside the cores. Ideas of Ref. Q were cast into a 
quantitative theoretical framework in Ref. where the 
vortex line was treated as a long ballistic quantum chan- 
nel. It was shown that the electronic thermal transport 
through a vortex can be described within the framework 
of the Landauer approach in a manner similar to that 
used for heat transport in normal conductors™**. The 



vortex heat conductance can be expressed in terms of 
the "heat conductance quanta" as a certain amount of 
heat carried by each quantized conducting modei^. This 
yields kl = kqNl, where the quantum kq = irT/(ih) is 
the universal heat conductance per conducting mode in 
a normal metal (we use the units with ks = 1), and Ni, 
is the effective number of conducting modes in the vor- 
tex (for both spin projections of electron). Quantization 
of the heat conductance in normal conductors has been 
observed in several experiments (see, e.g., Refs. [TTIIT^ |). 

To estimate the number of modes Al we summarize 
the main properties of the vortex core spectrum. The 
CdGM spectrum, e M , is a function of the quantized (half- 
integer) angular momentum fi = 6fcj_, where b is the im- 
pact parameter, k± = y/kp — fej, fop is the Fermi mo- 
mentum, and k z is the momentum projection on the vor- 
tex axis. As /i changes from — oo to +oo, varies from 
— Ao to +Ao crossing zero when /i changes its sign. Here 
A is the gap value far from the vortex axis. At small 
energies, |e| -C Ao, the spectrum is ~ — ^Ao/(fcj_£), 
where £ = ?jVf/Ao is the superconducting coherence 
length, and Vf = hkp/m is the Fermi velocity. The 
inter level spacing u>q ~ Ao/(&f£) taken at k z — de- 
termines the excitation minigap ojq/2. For temperatures 
ujq *C T -C Ao, the number of modes is thusS ~ T/ojq, 
and the heat conductance of one vortex becomes 

k l ~ T 2 /hu . (1) 

This estimate can also be obtained using the Sharvin 
conductance with the group velocity v g — de^/hdkz ~ 
e^/hkp rather than v z ~ vf as in a normal tube. Com- 
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paring Nj_, with the corresponding Sharvin's number of 
channels in the normal wire of radius £, Nc,h ~ (^f?) 2 ! 
one sees that N L /N Sh ~ v g /v F ~ (T/T c )(fc F ^)~ 1 < 1. 
At extremely low temperatures the minigap in the spec- 
trum suppresses the heat transport exponentially: n oc 
exp(— luq/2T) when T <C too- 

The qualitative estimate of Eq. applies to the heat 
conductance along the vortex core in the samples infinite 
in the lateral direction, but as we will see below does not 
generally work for samples finite in the lateral direction. 
Recent remarkable progress in experimental technique 
made it possible direct heat transport measurements in 
mesoscopic normal-superconducting heterostructuresi^ 
of sizes comparable to the superconducting coherence 
length, opening, in particular a new route for the analysis 
of quantum heat transport through individual vortices in 
mesoscopic superconductors whose exotic vortex states 
are in the focus of current research^. This calls for a 
comprehensive quantitative theory of the heat transport 
in mesoscopic superconductors. 

In this paper we develop a theory for an electronic heat 
conductance along the normal channels (vortex cores) 
in the superconductors of the arbitrary lateral size and 
shape, focusing on the specific features that appear in 
the mesoscopic cylinders with diameters comparable to 
the coherence length £. We show that in finite supercon- 
ductors, especially those with the a transverse size of the 
order of a few coherence lengths £, the estimate Eq. JTJ is 
no longer valid. The quantum mechanics of quasiparticle 
excitations in mesoscopic superconductors turns out to 
be very sensitive to the boundary effects. In particular, 
the interplay between the Andreev reflection at the core 
boundary and the normal scattering of quasiparticles at 
the sample surface results in an oscillatory behavior of the 
energy levels as functions of the Fermi momentum and 
the sample lateral sizei^. These conclusions have been 
recently confirmed also for a graphene channel in contact 
with a superconducting environmentpit For a vortex in 
a mesoscopic sample, the amplitudes of level oscillations 
can well exceed the CdGM interlevel spacing. In this 
case, the new quasiparticle modes propagating along the 
vortex core appear which leads to a dramatic increase 
in the effective number of the conducting channels, Al, 
as compared to that in Eq. Q). We present a theoreti- 
cal analysis of the single-particle transport mechanisms 
in the mixed state of mesoscopic superconductors tak- 
ing into account normal reflections of the quasiparticle 
waves at the superconductor boundaries. We concentrate 
on the low-temperature T <^ Aq regime for the ballistic 
samples with the mean free paths, both elastic and in- 
clastic, being much larger than the sample dimensions. 
We further investigate the effects of the surface rough- 
ness on the excitation spectrum and, accordingly, on the 
heat transport. We show that the significant increase in 
the number of conducting modes caused by the spectrum 
oscillations is a generic feature of a broad class of surface 
imperfections. 

The paper is organized as follows. In Section [H] we 



introduce the model and consider the basic equations de- 
scribing the quasiparticle quantum mechanics for a vor- 
tex line placed in a mesoscopic superconducting cylinder 
with the surface of an arbitrary shape. A Landauer-type 
analysis of the heat transport due to the quasiparticle 
modes propagating along the vortex core is carried out in 
the Section ITTT1 In this section we also present the results 
of numerical calculations of the vortex heat conductance 
in an ideal mesocopic cylinder. Then in the Section IIVI 
we investigate the effect of the surface roughness on the 
quasiparticle spectrum. 



II. QUASIPARTICLE SPECTRUM FOR A 
VORTEX LINE IN A MESOSCOPIC 
SUPERCONDUCTOR 

A. Semiclassical equations 

In this section we consider a vortex placed in a super- 
conducting mesoscopic cylinder with an arbitrary (x, y) 
cross section and derive an eikonal approximation of the 
Bogoliubov-de Gennes (BdG) equation that allows us to 
obtain the vortex-core states in such samples. Let the 
vortex line be parallel to the cylinder axis z. We choose 
the origin of the polar coordinate system (r, 8, z) of the 
(xy) plane at the vortex center. The surface of the sam- 
ple is then defined by the equation r = R(9). Since the 
momentum component along the vortex axis is a con- 
served quantity, we choose the particle- and hole-like 
parts of the wave function as (U, V) = (u, v) exp(ik z z) 
respectively. The BdG equations have the form 



Here 



H (p) A(x,y) 
A*(x,y) -JT *(p) 



Ho(p) = 7T- IP 
2m 



(2) 



E\ 



p = -ihV, Ex = E F -h 2 k 2 J(2m) = h 2 k 2 L /(2m), and E F 
is the Fermi energy. The order parameter profile A(x, y) 
is homogeneous in the z direction, 

A(x,y) = A 5 v (r)e^ . 

Here S v (r) is the normalized order parameter magnitude 
for a vortex centered at r = 0, such that 5 v (r) = 1 for 
r — > oo, e llf = (x + iy)/r is the vortex phase factor, and (p. 
is the part of the order parameter phase induced by the 
boundary effects; dp having no singularities inside the su- 
perconductor. For simplicity we start our analysis with 
the case of the weak external magnetic field and the ex- 
treme type-II superconductors where the vector potential 
A can be neglected. The effects of the magnetic field are 
discussed later in Section III CI 

Let us introduce a momentum representation: 
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d 2 pe lpr/h ^{p) , 



(3) 
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where p = p(cos9 p ,sm9 p ) — pp . The unit vector p 
parametrized by the angle 9 p defines the trajectory di- 
rection in the (a;, y) plane. Within the quasiclassical ap- 
proach the wave function in the momentum representa- 
tion assumes the form: 



*(p) = 



dse~ l(M - hk±)s/h ^(s,9 p ) 



(4) 



We look for solutions with |p| w hk±. For the slowly 
varying part of the wavefunction tp we obtain the quasi- 
classical (Andreev) equation along a trajectory specified 
by an orientation angle, 6 p , and an impact parameter, 
b, with respect to the vortex center. The distance along 
the trajectory is s = r cos(9 p — 9). The wavefunction for 
r = r(cos0, sin#) is derived from Eqs. ®: 



t ik^r cos(6 p —9) 



p- ti >^(rcos(9 p -9),9 p ) . (5) 



The wave functions should vanish at the sample surface: 



2tt 



^- j d9 p e ik ^ e)cosie "- e) ^(R(9)cos{9 p - 9),9 P ) = 



We introduce the angular-momentum expansion 



(6) 



,i IJ ,e p +ia z [e p +0(s,e p )]/2 
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where fi — —k±b = n + 1/2, and n is an integer. The 
function satisfies the equation: Hg^ = eg^, where 

d 

H = -ihV±a z — + AoS v (r)[& x s-a y b]/r 

+hk ± v s (s,b,9 p ) , (7) 

Vj_ = hk±/m, & x ,a y ,a z are the Pauli matrices, r = 
a/s 2 + b 2 , and 

h dip h ^ _ 
2m ds 2m 

is the projection on the momentum direction of the su- 
pcrfluid velocity component induced by the boundary ef- 
fects; we neglect this velocity assuming for simplicity that 
it is small as compared to the critical velocity at distances 
shorter than £ from the vortex axis. The function g^ is 
thus independent of 9 p : g^ = <7 M (s). This assumption 
does not change the main results of our work. 

The main contributions to the integral in Eq. © come 
from the vicinities of the stationary phase points, s'm(9 p — 
9) = /j,/k±R(9). We further assume that the r.h.s. of 
this equation is small, i.e. that the impact parameters 
for all trajectories contributing to the solution are much 
less than the system size. As a result we obtain two 



stationary points: 9 P ~ 9 + fi/k±R(9) and 9 P ~ ir + 
9 — fi/k±R(9). Furthermore, we neglect the corrections 
/i/k±R in the above expressions for stationary angles; 
this is valid provided (i) the additional phase shifts along 
trajectories caused by the finite impact parameters are 
small k±b 2 /R <C 1, and (ii) g^s) is a slow function on 
the atomic length scale. 

Near the stationary angle we can write cos(6> p — 9) ~ 
I- {B p - 9) 2 /2 and cos((9 p - B) ~ -1 + {B p - B - tt) 2 /2, 
respectively. The integral in Eq. © becomes 



d9 p e tk±R(e)cos ( "~ e) i>[R{9) cos(0 p - 9),9 P ] 
lk ^ R ^-™/ A ip{R{6),9) 



2tt 



Y k±R{9) 

+ e- tk±R(e)+m/4 ^{~R{9),9 + tt) 
Thus, the boundary condition Eq. © takes the form: 
4>{R{B),B) = - e - 2ik ± R W + ™/ 2 ip(-R{9), 9 + ir) . (8) 

To find the function <? M we can use the results of RefiS. 
Introducing 



- _ e *,(tf-if signb-farctanf)^ 



(9) 



we obtain: 



d 

ihV±a z —Wn + a x A 6 v Wu = ( e 



HV± b 

(10) 

The two linearly independent solutions for low energies 

Q 

:e*: 



p — D(s)/2gj arctan(s/b)ir z -|-i(T z [sign(b) — l]7r/4 



e D-\D\ +1 2 



V (l+7 2 )(e 2Z3 + 7 2 ) 
and W2fj,{s) = sign(fe) Wi„(— s). We denote 



e u + ija z \ A 



1 ' 



D(s) 



2A f , , w , 2A 
W[ j5 v {s)ds , A = — 



e~ D ^ds , (11) 



Here 



7(Ai) = ^vT [£(/i) ~ £] I e ~ D(s)ds - ( 12 ) 



2A 2 M [°° 6(s) DW iiu k F 
Hk ± V ± Aj ~ e dS ~ k ± 
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is the CdGM energy spectrum for an infinite sample. 
We take the function as a sum 

of two linearly independent solutions 

&!„ = e +*.(*f-*i»ign*-i««t«nf) {w Xll +W 2ll )/2 , 



G 



2fi 



3 +».(if -if signb-^arctanf) ^ _ ^ ^ 



The characteristic length scale of the function g^ is of the 
order of £. Assuming 7 <§; 1 these two solutions can be 
rewritten in the form: 



G 2M = e iS ^ D ^ 2 a z X . 



Substituting all the above expressions into the boundary 
condition and taking account of the fact that for small /1 
we have (p{R, 9) ~ (p{—R, 9 + it), we obtain 



iTT/i — ias 



A = 



where ag = 2k^R(9), Dg = D(R(9)). These expressions 
hold as long as e~ De <C 1. Let us now introduce the 
angular functions 
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>(*)=£■ 



, C, 



i- 1 e-w Cj (d)de . 



They satisfy 



;7Ci( 



- 2 7<* 



7Tj +C 2 ( 



— „-Dg+ia g 



d(0) 



7 is now considered as a function of \x which is trans- 
formed into an operator 



1 = 7(A) =7 -* 



Since c, 



2tt) 



-Cj(9) while ag and -D# are 2tt- 



periodic we obtain 



1 
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-) = o( 



, — Dg+iag 



,—Dg +7T -iag + 7 , 



ci(9) , 



ij Cl (9 + ir) = ( e - Ds+io,B + e- Dl >+"- la( >+T) Cl (9) . 

This set of equations demonstrates that the normal re- 
flection at the surface couples the two trajectories going 
in the opposite directions that are described by two in- 
dependent functions c\ (9) and c\ {9 + 7r) . Our derivation 
uses trajectories with small impact parameters and, thus, 



assumes that these equations are valid provided the an- 
gular harmonics c\ u vanish for \/j,\ > min{fcj_£, ^/k±R}. 

Following Ref. it is convenient to introduce a two- 
component function 



ci(6) 
d(9 + n) 



(14) 



defined for < 9 < it with the boundary condition c{9 + 
7r) = i<7 y c(9). The final equation for c reads 

hc = ec . (15) 

where the effective Hamiltonian is h = ho + V, 

h = e { -i^) ,V = -A-^J° (16) 



and 



vg = 1 (e 



-Dg-iag 



-D g+ „+i- 



(17) 



Equations l|15|) ~ (|17[l are the central result of this Sec- 
tion. They provide a semiclassical description of vortex 
core excitations taking into account the normal scatter- 
ing processes at the sample surface of a general shape. 
The normal scattering is introduced by the potential V 
through the factors e lae which result in mesoscopic oscil- 
lations of the spectrum as functions of momentum and 
sample sizei&. We employ these equations to study the 
effects of surface imperfections further in Section llVI 



B. Spectrum of an ideal cylinder 

For an ideal cylinder with a constant radius R{9) = R, 
the solution of equations (|15|) - (|17[l were obtained in 
Ref. [Ifj ; we reproduce them here for the sake of com- 
pleteness and convenience. The wavefunction is 



ifi6 



1 



(18) 



/2tt \e™» 

where jj, = n + 1/2, which gives the energy spectrum^ 
2e -D{R) 



A 



■ cos(2k ± R — 7r/i -f- 7r/2) , (19) 



Using a simple vortex core models 

5 v (r) = r/^Te v , (20) 
where £ v ~ £ is the size of the vortex core, one obtains 
£ M _ A* /C (a) 



A 



k±£ v K-i(a) 



afCi(a) 



2k ± R - 7r/i + -J .(21) 
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FIG. 1: The spectra of quasiparticles for a vortex in the 
cylinder with R/£ = 3.5 (a) and R/£ = 4.5 (b). The gap is 
approximated by Eq. I|20|l with £ v = £, /cf£ = 200. 



Here a = 2(£ v /t;)kp/k± and /C„ is the McDonald func- 
tion. Typical energy spectra of Eq. I|21|) for a vortex in a 
mesoscopic cylinder are shown in Fig.^ The minigap in 
the spectrum exists only for R > R c , where the critical 
radius is 



R< 




2fc F £ 



2fc F £ 



/Co (2&/0 



(22) 



_£o(2&/0 

The interlevel spacing wo for CdGM states is 

"°-fc^M2W0- (23) 

Taking for simplicity £ v = £ we get wq — 0.8143Ao/fc F £. 



Magnetic-field tuning of the spectrum 
oscillations 



In this section we consider the influence of an external 
magnetic field on the spectral characteristics and show 
that the quasiparticle spectrum can be tuned by changing 
the magnetic field H = Hzq applied to the sample. 

With the account of the magnetic field, the CdGM 
energy is renormalizedi^, — > + IMjJ c /2 where to c — 
eH/mc is the cyclotron frequency. However, the most 
important is the effect of the magnetic field on the gap 
profile. Far from the core it becomes 



n 2 eH 2 r 2 
20g 



(24) 



where 4>q — nhc/e is the magnetic flux quantum. We 
assume that the supercurrent density at the sample edge 
is less than the critical current so that the role of states 
bound to the edge^S can be neglected for low energies. 
The change in the gap profile results in an additional 
small renormalization of the minigap and, what is more 
important, it also affects the amplitude of level oscilla- 
tions determined by the factor exp(— D(R)). Indeed, we 



find from Eqs. (TTJ, HU 
D(R) 



T7 



2R 

T 



3R 



where h = H/H and H — (fro/nR 2 . The field H is of 
the order of the critical field Hi = (4>o/ttR 2 ) ln(i?/£) at 
which the vortex becomes energetically favorable in the 
sample. The ratio of the oscillation amplitudes taken, 
e.g., for h = 3 and h = 1 is given by the expression: 

e -D{R,h=i) - ex P {skxRj ' 

For k± = kp and R/£ = 4 this estimate gives us a fac- 
tor ~ 2, which shows that the amplitude of oscillations 
can be varied significantly by changing the magnetic field 
in a reasonable range. The critical radius R c at which 
zero energy modes appear is determined by the condi- 
tion D(R C , H) ~ ln(fc F £) and does also depend on H. To 
conclude here, one observes that the energy level oscilla- 
tions and the number of propagating quasiparticle modes 
depend strongly on the external magnetic field and thus 
can be effectively controlled in the experiments. 



III. THE LANDAUER APPROACH TO THE 
VORTEX HEAT CONDUCTANCE 

A. Semiclassical scheme 

Now we generalize the Landauer approach of Ref. 
developed for the electronic heat transport in supercon- 
ductors with infinite transverse dimensions to mesoscopic 
samples and derive the along-the- vortex heat current. A 
general expression for the energy current along z is 



h 



dk z 
■Km 



(hk z - ~ c ^z) v>nk x n(e,j,) 



(25) 



Particles with the distribution n(e) carry the energy +e 
while the holes with the distribution 1 — n(—e) carry 
the energy — e. If the electrodes are in equilibrium, 1 — 
n(— e) = n(e) at each electrode. The sum is taken over 
the states with e M > 0. Using the BdG equations one can 
derive the important identity for the wave functions of 
the bound states u il \ tz , v^u z belonging to the eigenvalues 
e M of a given momentum k z along the z-axis^ 



he 



d 2 r 



For the CdGM case this identity reveals a huge cancel- 
lation in its left hand side where each term is by the 
factor fc F £ larger than the right hand side. Within the 
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quasiclassical approximation the l.h.s. vanishes for an 
infinite vortex line as a direct consequence of an approxi- 
mate electron-hole symmetry. Note that such a cancella- 
tion does not occur in a finite-thickness slabS where the 
CdGM states are not truly localized. Making use of this 
identity reduces the expression for the heat current to 



u = 



Ef , \ den dk z 
J £ ^W z ^n 

V -fe F 



(26) 



Equation (|26[1 shows that, in a contrast to the electri- 
cal current, the energy flow is determined by the group 
velocity of excitations and this is why the conductance 
in Eq. Q is much smaller than Kg^ = kqN^, where 
kq = TrT/(Sh) is the heat conductance quantum defined 
in the Introduction. 

The quasiparticlcs in the injector lead (left lead) hav- 
ing energies coinciding with the core-state levels pene- 
trate into the core and are then transmitted through the 
vortex to the right lead with a probability P Mi / £i . In the 
same manner as it was done for the charge transport 
in mesoscopic conductorsiSiiiiSi we conclude that excita- 
tions with a positive group velocity v g = de^/Hdk z have 
the distribution n^P^k, where is the distribution they 
had in the left lead. Similarly, excitations with a negative 
group velocity possess the distribution urP^^z where tir 
is their distribution in the right lead. Since e^ is an even 
function of k 7 we obtain 



Defining the heat conductance as 



dk 2 



dk z 



(27) 



h = k(Tl - T R ) 



where Tl^r are the temperatures of the left and right 
electrodes respectively, one finds 



irTh 



E 



1 de,. 



P 



de, 



dk 2 



dk. 



(28) 



We assume that Tl — Tr, <C Tl,r, ~ T. Using the conduc- 
tance quantum kq, one can introduce the effective num- 
ber of transverse modes propagating along the vortex line 
N e f f — k/kq. Later in this section we neglect the scatter- 
ing of excitations at the interface between the leads and 
the superconductor assuming that they are transmitted 
freely through the vortex, i.e. with P^fc^ = 1, deferring 
the effects of scattering till Section IIII Bl 

For an infinite sample R 3> R c in the limit loq <C T <C 
Ao the sum in Eq. (|28|l can be replaced with the inte- 
gral over dfj, for the states with positive energies. Since 
dn/de = -{AT)- 1 cosh~ 2 (e/2T) we find 



K = 



9C(3) T 2 

27T HlUq 



Here £(n) is the Riemann zeta-function. Therefore 



W M (T) _«C(3)T 



2tt 2 



T 



which coincides with the estimate discussed in the intro- 
duction (see Eq. Q). At extremely low temperatures, 
T < uj : 



N<$(T) 



3 
4^2 



2TJ 



(29) 
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FIG. 2: The temperature dependence of the effective number 
of modes N e ff(T) calculated using Eq. with P M)fea = 1 
for two cylinders with R/£ = 3.5 (a) and R/£ = 4.5 (b). The 
gap is approximated by Eq. H20H with £„ — £, fci?£ = 200. 
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FIG. 3: The effective number of modes N e ff according to 
Eq. I|28|l with P^.fc^ = 1 as a function of the cylinder radius 
for various temperatures T/u>o = 0.3, 1.0, 3.0 (from bottom 
to top). The gap is approximated by Eq. I|20|l with = £, 
fc F £ = 200. 

In a mesoscopic cylinder with R < R c there appear 
new conducting modes due to the oscillatory behavior of 
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the energy levels. Indeed, since it is the absolute value 
of the group velocity that enters Eq. l|2^|) , both positive 
and negative slopes of the oscillating energy e^fcz) as a 
function of k z in Eq. H19|) or in Fig. ^ contribute equally 
to the heat flux. This contribution is finite even for very 
low temperatures and is determined by the number of 
crossings of the constant energy line e ~ T by the oscil- 
lating spectrum. Thus this number represents the num- 
ber of conducting modes N^l ~ N^N^, where N^ z is 
the number of crossings, as a function of k z , of a spectral 
branch with a given angular momentum /i, and is the 
number of the angular momentum branches that cross 
this line of energy e ~ T. For large number of modes, 
i.e., when the oscillation amplitude is larger than lvq we 
havei&iV fez ~ k F R(l-R/R c ) andiV^ ~ (2A /ujo)e~ 2R ^ . 
Finally, we get 



N 



{ f f ~ k F R(l - R/R c )(2A /cJ )e- 



-■2R/t 



(30) 



If R ~ £, then N ( e fj ~ N sh ~ (k F 2 whi ch coincides 
with the estimate of the number of transport modes for 
a normal metal cylinder. We see therefore that the os- 
cillations of the vortex-core energy spectrum in a meso- 
scopic cylinder lead to a drastic increase in the number of 
the conducting modes and, accordingly, to an enhanced 
electron thermal conductance through the vortex cores 
as compared to that in an infinite sample. 

The temperature dependence of the effective number 
of modes calculated using Eq. (|28|l is shown in Fig. |5| 
The effective number of modes for T 3> loq is 



N$(T)~N* 



27((3) T 



2tt 2 



LJ 



(31) 



Here N* is constant as function of temperature; it is 
small for R > R c . However TV* grows and approaches 

the value N { e f f in Eq. (|30[l if the cylinder radius decreases 
below R c and the amplitude of oscillation increases. This 
asymptotic temperature dependence is observed already 
for T > 0.5 wo- The initial steep increase in N e ff in 
Fig. [21 is associated with the crossover into the regime 
of large number of modes when temperature is increased 
up to T ~ loq. The constant value N* leads to a linear 
temperature dependence of the heat conductance k(T) 
in the range T > ujq. The linear in T contribution can 
exceed the quadratic term given by Eq. or by the 
second term in Eq. (|31|l . 

The effective number of modes calculated using 
Eq. (|28|l is shown in Fig. [31 as a function of the cylin- 
der radius. Its non-monotonic dependence on R agrees 
reasonably well with our simple estimate in Eq. H30[) . 



B. Exact quantum- mechanical analysis 

A disadvantage of the above semiclassical picture is 
that it does not provide a systematic way how to cal- 
culate the probability for excitations to penetrate into 



the vortex core modes, which is accompanied by reflec- 
tion of excitations from the vortex ends at the boundaries 
between normal reservoirs and the mesoscopic supercon- 
ductor. In this Section we derive the general quantum- 
mechanical expression for the heat conductance in super- 
conductors directly from Eq. (|25|l and demonstrate that 
it has the form of the Landauer formula with the proper 
account of the scattering processes. Next we find the 
exact solution of the BdG equations and calculate the 
vortex heat conductance numerically. 

We consider a cylindrical superconducting slab of the 
thickness d (Fig. 0J> confined between the two normal 
leads (reservoirs) . The normal-superconducting interface 
is assumed ideally transparent. The axis of a trapped 
vortex line coincides with the superconducting cylinder 
axis. To calculate the heat conductance we solve the 
scattering problem for electron-like (U) and hole- like (V) 
quasiparticle wave functions within the BdG theory. It 



r, 


R 


heat flow 
U 

N 






-> vortex core 




N 


s 



FIG. 4: Heat transport through a vortex line. 

is convenient to consider four types of linearly indepen- 
dent states: (i) the state (Ui,Vi) with an electron wave 
incident on a superconductor from the left lead in the 
absence of other incident electrons or holes; (ii) the state 
(U2,V2) with a hole wave incident on a superconductor 
from the left lead; (iii) the state (C/3, V3) with an electron 
wave incident on a superconductor from the right lead; 
and (iv) the state (U^, V4) with a hole wave incident on 
a superconductor from the right lead. 

In the left lead, the first (electron) state has the form: 



,i( P +i/2)e 



U, = 



/2tt 
e i(n+i/2)e 



Vi 



+ - 

e i( P -i/2)e 



^ Tuti, ss' Quips' ix)^ zu , (32) 



r U v,ss'9v^s'{r)e +lk "' z . (33) 



Here g u ,fj,s(r) and gv^sir) are the transverse electron and 
hole modes of a normal metal wire with an angular mo- 
mentum /i and a radial quantum number s respectively. 
The transverse modes for each /i satisfy the boundary 
conditions g u .,_ ls (R) — g v ,^s(R) — 0; they are normalized, 



iu,psy u ,fj,s' 



■ dr = 



Jv,fisg v ^ s ' r 



dr = 5* 
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and have following form: 

V2J |M+1/2 , (# +1/2|) r/i? 



fJa. 



fJS 



RJ\ft+l/2\+l [Js 
Tot / 1/2|) 



9v„ 



(|M+1/2|) 

r/R 



RJ\, 



-1/2| + 1 IJs 



:(|M-1/2|) 



(34) 
(35) 



Here J n is the Bessel function of n-th order, ji^ is the 
s-th root of the equation J„(ji™' ) ) = 0. The solution Eqs. 
(|32|l . (|33(l corresponds to the energy 



2m 



(ii^ +1/2|) /^) 2 - 



(36) 



The wave numbers of the scattered electron and hole 
waves are k zuv = k' zuv (fi, s' , k z ), where 

k' zu J^ S ',k z ) = yfc|-(jl! A1±1/2|) /i?) 2 ±2m £ /^ . 

The factors r UUtSS i and r uv>ss i are the matrices of nor- 
mal and Andreev reflection amplitudes for the processes 
with large and small transfers of longitudinal momen- 
tum, respectively. They are non-diagonal in s and s' 
which accounts for scattering with a small change in the 
transverse momentum caused by an inhomogeneous vor- 
tex order parameter profile in the radial direction. The 
change in the radial momentum at a given energy e leads 
to a change in the longitudinal momentum such that the 
scattered waves have momenta k zu and k zv slightly dif- 
ferent from the momentum k z in the incident wave. 

Similarly, the state with an incident hole wave is 



2tt 

i{n-i/2)e 
V 2 = ^—Q, Kll ,(r)e-^ z 



(37) 



/2tt 
e i(fi-i/2)e 



The energy of the incident hole is 

h 2 r 



9v,fis \ 

Y. r ™,ss>gv,»Ar)e +lk '^ . (38) 



e(lJL,s,k z ) = - 



2m 



(39) 



The transmitted state (iii) in the left lead is 

e Hp + l/2)0 



'2tt 



s' 



(40) 



with the replacements t U u,ss' ~~ > t VUtSS 
tw,ss' • The factors t uu ss * and tw,ss 



The waves U4, V4 in the state (iv) have the same form 

and t U y^ ss ' ► 
/ describe the 
particle-particle and hole-hole transmission processes 
with a small change in the longitudinal momentum while 
the factors t uv _ ss i and t vu ^ ss i describe the transmission 
with a large transfer of the longitudinal momentum in- 
volving normal scattering processes. 

BdG equations (0 conserve the quasiparticle current 



1 

2m 



u*(-ihV - -A)u - v*(-ihV + -A)v + c.c. 
c c 



For the incident and transmitted waves Eqs. (|32|) . 13311 . 

etc., this conservation law yields the sum rules 



Tuu,ss' I &zv \^uv,ss' I 



~^~fezu\tuu,8S f I ~l~ k zv \t UVj88 t I ) — k z , (41) 

{^zv\ r vv,ss' I + k zu \r VU]SS ' I 



(42) 



In the quasiclassical approximation, the transverse mo- 
mentum is conserved, therefore 

At the same time, the processes with large longitudinal 
momentum transfer arc prohibited 



' uu.ss' — ' vv,ss' 



— tuv,ss f — tvu,s 







However, in the exact BdG formalism used in our nu- 
merical calculations, we keep the full set of probabilities 
including t uu ss / 7 t vv ss i 7 t uv ss ' 7 and tvu,ss'- 

The heat flow through the cylinder cross section in the 
left lead can be calculated summing up the states (i)-(iv): 



Is = J2 [ rdrd9 ^{ ^■^{UiyU 1 }n L (e)~e-Im{V 1 *VV 1 }-(l-n L (-e))] 

fl.S 

+ [e ■ lm{U* 2 VU 2 }n L {s) - e ■ lm{V 2 *VV 2 } • (1 - n L {-e))) 
+ [e ■ lm{U;WU 3 }n R (e) - e • Im^Wa} ■ (1 - n R {-e))} 

+ [e-lm{U:\7U i }n R (e)-e-lm{V:VV i }- {l-n R (-e))} } . (43) 
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The sum is taken over the propagating quasiparticle 
states with Imfc 2 = and e > (we thus take into ac- 
count both incident electrons and incident holes). 

Using the sum rules Eqs. (|41ll . ((42(1 and transforming 
to the energy integral, the heat current becomes 
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irh 



e[n L {e)-n R (e)]T(e)de, (44) 



where the transfer factor 

El , 2 2 1 2 2 

\^uu,ss' I H~ \tuv,ss' I ~t~ I^-du.ss' I H~ |^v^,ss ; I 

. (45) 

is the sum of partial transmission probabilities. Finally, 
the dimensionless heat conductance (the effective number 
of modes) becomes 



Neff = 



de 



(46) 



Equations 1(44(1 . 1)45(1 . 1(46|) are the exact quantum- 
mechanical expressions for the heat conductance. They 
correspond to the Landauer formulaic in which the trans- 
mission probabilities are defined through the full scatter- 
ing problem for excitations in superconductors. 

The semiclassical approach developed in the previous 
section assumes that each particle or hole, incident on 
one of the superconductor/normal- metal interfaces, with 
the energy and momentum corresponding to the one or 
another propagating vortex core mode is subsequently 
transmitted along this channel with a probability Pi to 
the second normal lead: 



provided the direction of the quasiparticle momentum 
along z coincides with the direction of the group velocity 
in the vortex core, and 



in case where the direction of the quasiparticle momen- 
tum along z is opposite to the direction of the group ve- 
locity. The index i denotes the set of quantum numbers 
H,k z ; the sum is taken over such quantum numbers sf 
and s~ which satisfy the energy conservation equations 
e n[k z {sf )] = e for increasing or decreasing branches of 
e M as functions of k 



respectively (we emphasize that 
the quantum number s is regarded as continuous in the 



semiclassical approximation) . 
can write Eq. 1)45(1 as 



Using t v 



one 



dk z 



6(e - e t _ l (k z ))dk z 



As a result, Eq. (|4^|l coincides with Eq. (|2"7}> . Under 
the assumption of complete transmission P^k z = 1 the 
transfer factor T(e) is simply equal to the number of 
quasiparticle modes propagating along the vortex core 
for a given energy e. 




0.5 k' /k c 

zu F 



0.5 k' /k c 1 

zv F 



FIG. 5: (a), (b) The transmission Ituul 2 , |£ ut ,| 2 and (c), (d) 
reflection |r uu | 2 , |r uu | 2 probabilities, respectively, as functions 
of momenta k' zu and k' zv of scattered waves for the incident 
electron-like quasiparticle wave with fi = —1/2, k z /kF = 0.9, 
and £ = 0.49ljo- Parameters are: R = 4.5£, d/£ = 16. 



C. Numerical results 

Here we present the results of numerical calculations 
of the vortex heat conductance for a typical parameter 
A /E F = 0.01 (k F £ = 200). The vortex core was ap- 
proximated by Eq. 1(2011 with = £. The critical radius 
estimated by Eq. ((22(1 for the chosen set of parameters is 
R c ~ 3.96 £. The calculations were carried out for two 
values of R — 3.5 £ < R c and R — 4.5 £ > R c , and for sev- 
eral values of superconductor lengths d/£ varied from 16 
up to 2048. The corresponding energy spectra for prop- 
agating quasiparticle waves in a cylinder of an infinite 
length tabulated using Eq. l(2T)) are shown in Fig. ^ 

We calculate the scattering matrix of the normal 
modes incident on the mesoscopic superconductor with a 
trapped vortex line using a representation of the BdG op- 
erator in the truncated basis of the normal-metal eigen- 
functions. We solve the system of linear equations which 
correspond to the boundary conditions at the normal- 
metal/superconductor interfaces and to the conditions 
at the infinity (radiation condition). Typical scattering 
matrix coefficients as functions of the momentum of scat- 
tered wave are plotted in Fig. [SJ The coefficients display 
peaks around k zu ^ v = k z . The energy dependence of the 
transfer factor T(e) calculated using Eq. ((45(1 for two 
lengths d of the vortex line is shown in Figs. The lin- 
ear in e/u>o term in the transfer factor corresponds to the 
increasing number of CdGM modes with \fi\ujQ < e which 
contribute to the heat transport. 

The temperature dependence of the heat conductance 
calculated using Eq. 1(461) and plotted in Fig.[7|is qualita- 
tively similar to that obtained on the basis of the semi- 
classical model (see Fig. |5J). In particular, the slope of 
the linear asymptotic dependence (T > ujq) of N e ff vs 
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FIG. 7: The effective number of modes iV e //(T) calculated 
FIG. 6: (Color online) The transfer factor T(e) for different numerically as a function of temperature for the vortex line 
vortex line lengths d (blue line is for d/£ = 16 and red line length d/£o = 16 in cylinders with R/io = 3.5 (a) and ii/fo = 
is for d/£ = 2048) in the cylinder with R/£ = 3.5 (a) and 4.5 (b). 
R/t = 4.5 (b). 



temperature is indeed well described by the analytical 
approach. However, the value N e ff calculated numeri- 
cally is about seven times smaller than that evaluated 
within the semiclassical approach with P ll ,k z = 1- Such 
deviation is obviously caused by a finite probability for 
excitations to penetrate into the vortex core states at 
the normal- lead/superconductor interfaces. The conduc- 
tance peak at T* ~ 0.2wo which is most pronounced for 
the radius R > R c (see Fig. \J}p) is associated with the 
contribution of the lowest oscillating CdGM level to the 
transfer factor T(e). 

In our numerical calculations we do not observe the 
power-law dependence of the thermal conductance on the 
thickness d predicted in Ref. j9j| which is associated with 
the ballistic quasiparticle trajectories passing close to the 
vortex axis. This power-law behavior is expected to dom- 
inate only for much shorter d; the transport associated 
with such trajectories appears to be closed already for 
our set of material parameters and d > 16£ since it would 
result in N eff < 10~ 3 . 



IV. VORTEX IN A CYLINDER WITH A 
ROUGH SURFACE 

In the previous sections we calculated the electronic 
spectrum and the heat transport for a sample with a per- 
fect surface. One can expect, however, that surface im- 
perfections and roughness would reduce the amplitude of 
spectrum oscillations. Making use of the approach devel- 
oped in Section[n]we analyze here how surface roughness 



affects the vortex-core spectrum and the heat transport 
mediated by the vortex core states. We show that there 
exists a broad class of surface imperfections which do 
not destroy mesoscopic spectrum oscillations completely; 
in this case the heat transport through vortices still re- 
mains significantly enhanced as compared to that in the 
bulk system. 

Let us treat the potential V in Eq. (|16fl as a perturba- 
tion assuming the cylinder radius to be sufficiently large. 
The solution to this equation within the zero-order in 
V is a set of orthogonal and normalized eigenfunctions 
Eq. (|18fl . The first-order correction to the energy levels 
Se = e — e(/z) is 

A f 

Se = -4- / e- De cos \a e - it fi + Ti/ 2} d6 . (47) 
ttAJ 



Consider several examples. Assume first that the 
cylinder has a smooth surface R(9). The main con- 
tribution to the integral in Eq. (|47|) comes then from 
the stationary phase angles 0j defined by the condition 
da g /d6\ e . = 2k ± dR/d6\ e . = 0. In the limit k ± R" > 1 
we thus obtain 



a ^-y Trfexiimoi 

x cos \ ae 3 - 7T/x + n/2 + — sign (R (6j)) . 

For a simple realization R(0) = R + acos(m9) 1 where 
a <C £, the correction to the energy takes the form: 

Se = (2A /A) e~ D ° cos(a - irfi + ir /2)J (2k ± a) 
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where ao = 2k±Ro and Dq 
1 we obtain 

Se = (2A /A)(irk ± a)- 1/2 e- Do 

x cos(ao — 7r/i + tt/2) cos(2k±a — 7r/4) . 

The amplitude of oscillations is decreased by a factor 
{kpa)^ 1 / 2 <C 1. Nevertheless, the amplitude Se ^> loq as 

long as e' D ° > (fc F a) 1/2 (fcf Note that the r - h - s - 01 
this inequality is a small parameter. 

Let us now assume that the cylinder has a rough sur- 
face such that the radius is R(9) = R a + a(9), where a(9) 
is a random function of sample realizations with a zero 
ensemble average (a) = having a Gaussian distribution 



D(Rq). In the limit k±a ^ This average vanishes for kp<J ^> 1. 



However, the energy level fluctuations are consider- 
able. To prove this we calculate the mean square en- 
ergy fluctuation assuming a Gaussian distribution for 
ai = a(9i) and a 2 = 0(^2): 



102(01,02) 



1 



2vrcrVl - K 2 



exp 



2(1 - K 2 )a 2 



wi(a) 



27r<7 



exp 



a 

2^2 



Here a <C £ is the mean square value: (a 2 ) = a 2 . Aver- 
aging of Eq. (|T7|) gives 

(Se) = ^-e- D "cos(a -irn + ir/2) I ' * \e 2ik ^)d6 

7rA ,/n 



2A 



A 



-e D ° cos(ao — 7r/i + 7r/2)e 



-2k, 



where r = \/a\+ a| — 2Ka\ai and AT = A(#i — fe) = 
(a\ai)ja 2 is the correlation coefficient. This correlation 
coefficient should be an even and 2-7r-periodic function of 
angle: K{9) = K{-9) = K{6 + 2vr). For 6 1 = 9 2 we get 
K = 1. The correlation coefficient decays as we increase 
\9-\ — 9o\ value from to tt. We obtain 



(Se 2 



A 2 

_i±o_ e -2r> 



7T 2 A 2 



(cos(ae 1 — 71-jU + 7r/2) cos(ag 2 — 7T/i + 7r/2)) d9\d9 2 



A 2 

^0 _ e -2D 





27T 27T 



2tt 2 A 2 



cos 



(2ao) / e 2lfc -L( a i+ a 2)\ -I- / e 2ifcj_(ai-a 2 ) 



d6»i^2 







For k±a 1 the second term is only important. Using 
the characteristic function for a Gaussian process 



gi(«l 01+11202) \ _ e -^(T 2 (tl?+«2+2-R"'"lM2) 



where wi = — U2 = 2/cj_ we get 

2tt 27T 



Expanding the correlation coefficient for small angles, 
K(9) ~ 1 - |A^(0)|6» 2 /2, and introducing the correla- 
tion radius £ = R o /y/\K"{0)\ we find 



(fe 2 



A 2 ^ 



A 2 k ± (tR (27t) 1 / 



e -2D 



(48) 



This expression is valid if the range of angles contributing 
to the integral is small, 9 — l/akj_^\K"(0) \ < 1. Taking 



t ~ cr we obtain 



A,, 



--Do 



(49) 



We see that the mean square fluctuation is much larger 
than the CdGM interlevel spacing if the sample radius is 
not very large, e~ D ° > (k F R ) 1 / 2 (k F ^)^ 1 . Note again 
that the r.h.s. here is a small parameter. 

When the correlation radius increases and sample di- 
mensions are small i?o < £/k±a we can take K ~ 1 for 
all 9. As a result 



^JSe^) ~ (A 72/A)e 



--Dn 



(50) 



In this limit the fluctuations of energy levels are universal 
in the sense that they do not depend on I and cr, though 
do depend on the sample size. 

To summarize, we see that imperfections on sample 
surface do not always remove the spectrum oscillations. 
We have demonstrated that the number of propagating 
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quasiparticle modes responsible for the heat transport 
in a mesoscopic sample can still be essentially large for 
smooth surfaces with comparatively big variations in the 
radius and also for Gaussian imperfections with fluctua- 
tions as large as a kp 1 . 

V. CONCLUSIONS 

We have investigated the electronic thermal transport 
along the flux line traversing a mesoscopic superconduct- 
ing sample having the shape of a cylinder with the ar- 
bitrary cross section. We have shown that the normal 
quasiparticle scattering at the sample boundaries results 
in the essential increase in the number of modes propa- 
gating along the vortex core. We have found that the 
semiclassical version of the Landauer theory provides 
an adequate qualitative description of the heat trans- 
port carried by these modes. We have verified this con- 
clusion by the detailed numerical analysis of the exact 
quantum-mechanical scattering problem described by the 
Bogoliubov-de Gennes theory. For small-radius cylin- 
ders, such that R < R c , the minigap in the vortex spec- 
trum is suppressed and the conductance becomes a linear 



function of temperature in the range T > luq. This linear 
term in the vortex conductance increases with the de- 
crease in the cylinder radius; it can strongly exceed the 
quadratic term Eq. which usually dominates in the 
bulk system at u> Q < T < T c . 

We have analyzed how imperfections at the cylinder 
surface affect the quasiparticle spectrum. We have found 
that there exists a wide class of surface imperfections 
that give rise to only partial reduction in the number of 
the conducting modes mediating the heat conductance, 
even if the dimensions of imperfections exceed the atomic 
scale kp . 
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